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ABSTRACT 

We have determined the relativistic hght deflection of the quasar J0842+1835 as Jupiter 
passed within 3.7' on 2002 September 8, by measuring the time delay using the VLBA and 
EfFelsberg radio telescopes at 8.4 GHz. At closest approach, General Relativity (GR) predicts a 
radial (static) deflection of 1190 /larcsec, and a tangential (retarded) deflection in the direction of 
Jupiter's motion of 51 /iarcsec. Our experiment achieved an rms position error of < 10 /iarcsec, 
and measured this retarded deflection to be 0.98 ± 0.19 (rms error) times that predicted by GR. 
The increased positional accuracy for this VLBI phase referencing experiment was achieved by 
using two calibrator sources. Comments on the interpretation of this experiment, are given. 
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1. The Gravitational Deflection of Radio 
Waves by Jupiter 

Einstein solved the equation of light propaga- 
tion in the field of a static body, and predicted the 
deflection of light at the limb of the Sun of 1.75" 
(Einstein 1916), and this relativistic deflection was 
measured for the flrst time in 1919 (Dyson, Ed- 
dington & Davidson 1920). The deflection pre- 
diction by Einstein is equivalent to measuring the 
time delay of light traveling in a static gravita- 
tion field (Shapiro 1964), and was confirmed in a 
series of experiments, the most precise of which 
was made with Very Long Baseline Interferom- 
etry (VLBI) (Lebach et al. 1995; Eubanks et 
al. 1997). Kopeikin (2001, 2003) generalized the 
problem of light propagation in the gravitational 
field of arbitrary moving bodies, and showed that 
the Lorentz-invariant expression for the relativis- 
tic time delay of light (and light deflection), to all 
orders in v/c, depends on the retarded positions 



of the moving bodies as defined by the retarded 
Lienard-Wiechert solution of the Einstein gravi- 
tational field equations. For Jupiter, as an ex- 
ample, the Lorentz-invariant relativistic time de- 
lay between two telescopes (#1 and #2), from a 
source of light at infinity (quasar), is given by the 
equation (Kopeikin 2001) 



A = 



2GMj 



K-V,;(si) 



In 



rijjsi) + K-ri,7(si) 
r2j(s2) + K-r2j(s2) 



(1) 



where G is the gravitational constant, c is the 
speed of light, Mj is the mass of Jupiter, K is the 
unit vector from the barycenter of the solar system 
to the quasar, and vj(si) is the coordinate veloc- 
ity of Jupiter taken at the retarded time Si. The 
retarded times Si and S2 are calculated accord- 
ing to the gravity null- cone equations, as defined 
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by the Lienard-Wicchert retarded solution of the 
Einstein gravity field equations which contain a 
constant in the wave operator that we denote as 
Cg (Kopeikin 2003; Will 2003). The retarded times 
are 

si = ii-— |xi(ti)-xj(si)| , (2) 

S2 = i2-— |x2(t2)-Xj(s2)| , (3) 

where ti, t2 are times of arrival of the radio sig- 
nal from the quasar to the first and second VLBI 
stations respectively, having coordinates xi(ti) 
and X2(t2), and xj(s) is coordinate of Jupiter 
taken at the retarded time. The distances be- 
tween each telescope and Jupiter arc given by 
rij = \rij\, r2j = \r2j\, rij(si) = Xi{ti) -xj{si) 
and r2j(s2) = X2(f2) - xj(s2) , all of which de- 
pend implicitly on the value of Cg. GR predicts 
that the constant, Cg, equals the speed of light 
c. If Cg ^ c, then Eqs. (1) to (3) would not be 
invariant with respect to the Lorentz transforma- 
tion. Hence, the experiment measures the numer- 
ical value of Cg as a test of the Lorentz invariance 
of the Einstein equations, and is an indirect mea- 
surement of the speed of propagation of gravity^. 
Standard data analysis of VLBI observations uses 
the barycentric time t and the barycentric space 
coordinates x = {x,y,z). Hence, all quantities 
and variables in Eqs. (1) (3), and hereafter, are 
expressed specifically in these coordinates. For 
example, vj = dx/dt is the orbital velocity of 
Jupiter with respect to the barycenter of the solar 
system, and vector K defines the direction to the 
quasar as measured by an observer at rest with re- 
spect to the barycenter (see Kopeikin (2001) and 
Kopeikin & Fomalont (2002) for more detail). 

Relativistic time delay of light in the gravita- 
tional field of moving Jupiter was measured in 
a VLBI experiment on September 8, 2002 when 
Jupiter passed within 3.7' (^ 14 jovian radii) of 
the bright radio quasar J0842+1835. If Jupiter is 
at a small angular distance 6 from the quasar, Eq. 
(1) can be expanded in a Taylor series around the 

^Thc radio wave from the quasar docs not propagate along 
the gravity null cones (Eqs. (2) and (3)) and, therefore, 
Cg can not be physically identified with the speed of light 
(Kopeikin 2001, 2003), as was done, for example by Asada 
(2002) 



time of arrival^. Then, the linearized relativistic 
time delay, A{Ti^rn) between two radio telescopes 
from a VLBI network having numbers I and m and 
separated by a distance 6;^^) can be expressed to 
sufficient accuracy as 

A(t;,„) = 3 

A CgO ) e ^ CgO^ \ ^> 

where vtj = K x (vj x K) is Jupiter's velocity 

in the plane of the sky, n is the tmit vector from 
Jupiter to the radio source in the plane of the sky, 
and rmj is the distances between Jupiter to the 
mth telescope. The angle 9 is defined by the re- 
lationship (Kopeikin (2001), see Eq. (9)), cos = 
-K • p, where p = (xi - Xj(ti))/|xi - Xj(ti)|, 
and is calculated from the accurate JPL ephemeris 
(Standish 2000). The first term in the large brack- 
ets is the radial deflection term^, and the second 
term in brackets is the tangential deflection which 
is in the direction of the velocity of Jupiter with re- 
spect to the barycenter of the solar system, vtj, 
in the plane of the sky. The velocity-dependent 
terms in Eq. (4) result from the Taylor expansion 
around the time of arrival of ti of the retarded 
arguments in Eq. (1). In the analysis of the ex- 
perimental data, we used Eq. (l)-(3), rather than 
the approximate form in Eq. (4). This misunder- 
standing of the analysis process has led to the er- 
roneous claim (Samuel 2003) that no v/c terms 
are observed with this deflection experiment (see 
Kopeikin (2003) for more detail). 

The constant Cg, shown explicitly in Eq. (4), is 
associated in our model of the experiment with a 
speed of gravity ( scalar and vector modes) propa- 
gating along the gravity null-cone Eqs. (2) and (3) 
in accordance with the physical interpretation of 
the Lienard-Wiechert solutions of Einstein's equa- 
tions (Kopeikin 2001, 2003)^^. We introduce a fit- 

^In the general case when angle 6 is not small the accurate 
retarded Eqs. (l)-(3) must be used. 

^This equation treats more precisely the radial deflection 
than the corresponding equation from (Kopeikin 2001). 

^The tensor modes (free gravitational waves) also propagate 
along the gravity null cones but they can not be detected in 
the neaj" zone of the solar system. Hence, the experiment 
can not be viewed as a detection of gravitational waves 
which are too faint to be seen with present VLBI accuracy 
(Kopeikin et al. 1999). 
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ting parameter S ~ c/cg — 1 to measure the dif- 
ference between the two speeds. For GR, 5 = 0. 
Other theoretical models of this experiment (Will 
2003; Asada 2002) predict the same time delay 
given in Eqs. (1) and (4), but only to the first or- 
der of v/c. However, these different formulations 
lead to a different interpretation of Cg . See §4b for 
further discussion. 

The order of magnitude of the deflection pre- 
diction on September 8 for the 6000-km telescope 
separation is a delay of 115 psec (deflection of 
1190 /xarcsec) for the radial term, and a delay 
of 4.8 psec (deflection of 51 /larcsec) for the re- 
tarded term at the point of closest approach. Al- 
though f j/c « 4 X 10~^, the additional factor of 
1/6 amplifies the retarded term so that it is 4% 
of the radial term^. A previous close passage oc- 
curred in 1988 March 21 (Treuhaft & Lowe 1991) 
and the radial deflection term was measured to 
an accuracy « 15% in accordance with GR. With 
improvements over the years in VLBI techniques, 
sub-milliarcsecond positional accuracy were now 
attainable, and measurement of the retarded term 
was feasible. 

2. The Experimental Strategy 

The VLBI experiment to measure the deflection 
of light of the quasar J0842-I-1835 by Jupiter was 
conducted during five days, centered on Septem- 
ber 8, 2002. The radio array consisted of eleven 
telescopes: ten 25-m diameter telescopes of the 
Very Long Baseline Array (VLB A), plus the 100- 
m diameter telescope at Effelsberg, Germany. The 
observing frequency was 8.45 GHz, and ten hours 
were used for observations on each of 2002 Septem- 
ber 4, 7, 8, 9, 12. 

An array measures the difference in arrival 
time, the delay, for a quasar signal to reach each 
of the telescopes. Using a widely spaced array, 
extremely high positional accuracy can be ob- 
tained. However, this delay is contaminated by 
many effects, both internal to and external to 
each telescope system. To remove these effects, 
radio astronomers use the technique of phase ref- 
erencing (alternating observations of two sources) 
whereby observations of one source (the calibra- 
tor) are used to determine the delay errors asso- 

'The amplification is due to the time lag between present 
and retarded positions of Jupiter (Kopeikin 2001). 



dated with the target source (Beasley & Conway 
1995). If the switching time and source separa- 
tion between calibrator and target are less than 
the temporal/spatial variation scale sizes, the rel- 
ative position between the two sources can be ac- 
curately measured. However, the propagation de- 
lay along the quasar paths through the ionosphere 
and troposphere may be so variable in time and 
angle, that even fast switching will not completely 
remove these propagation changes, and relative 
position accuracies < 30 /xarcsec are difficult to 
obtain even with VLBI techniques. 

After several sets of test observations in early 
2002 (Fomalont & Kopeikin 2002), we decided 
to sequentially observe J0842-I-1835 with two 
known nearby calibrators: J0839-f 1802, separated 
0.82° in position angle -132°; and J0854-h2006 
(=OJ287), separated by 3.36° in position angle 
63°. The sources are identified as quasars with 
high redshift, and they lie on a nearly straight 
line, which was exploited in the calibration proce- 
dure. The observing sequence was identical for the 
5 observing days: We cycled observations through 
the three sources in 5.5 min, with about 100 such 
cycles per day, a reasonable compromise between 
removing the temporal effects with obtaining good 
accuracy of the measured delay. This observa- 
tional technique was an improved variant of that 
used for the solar bending experiments in 1974 
and 1975 (Fomalont & Sramek 1976). The jovian 
magnetosphere also produced a radial deflection of 
the radio waves (directed inward to Jupiter) and 
this anomalous bending was of concern. However, 
calculations suggested that this bending would 
be signiflcantly less than the retarded deflection 
term, and special observing techniques (observ- 
ing at two frequencies simultaneously) would have 
reduced the sensitivity of the experiment to the 
gravitational bending. The magnetosphere bend- 
ing is discussed in more detail below. 

3. The Data Reduction 

3.1. The Radio Interferometer Response 

A radio interferometer measures the complex 
spatial coherence function of the electromagnetic 
radiation fleld C/.m e*^'^'^' ™ at a frequency v be- 
tween two telescopes denoted by I and m. If 
the intercepted electromagnetic radiation fleld is 
dominated by a small-diameter radio source, then 
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the response of the interferometer, denoted as the 
complex visibihty function Vi^: 



IS 



Vi, 



GiG* e^^^^'-^^-) 



(5) 



The complex gains of the telescopes, Gi, contain 
the amplification and phase shifts that are intro- 
duced by each telescope system (* denotes com- 
plex conjugate), and the residual time delay (ob- 
served - model) of the signal from the radio source 
to each telescope is denoted by At; . For an array 
of A'' telescopes, the visibility function is sampled 
simultaneously with N{N — l)/2 interferometers, 
55 pairs for an eleven-element array. 

The accurate delay model was calculated using 
the Goddard Space Flight Center (GSFC) CALC 
software package version 9.1^, which uses the most 
recent parameters associated with the earth ro- 
tation and orientation, nutation, the terrestrial 
reference frame, and the radio source positions 
defined on the celestial reference frame (Ma et 
al. 1998). We have slightly modified the CALC 
package to incorporate the gravitational deflection 
from all significant solar system bodies including 
the retardation term, as given by Eq. (1). We have 
assumed that the PPN parameter 7 = 1, because 
it has been measured to 0.1% accuracy (Lebach 
et al. 1995; Eubanks et al. 1997). We could not 
improve this value in our experiment, and its un- 
certainty produces a deflection offset < 1 /larcsec, 
much less than the experimental accuracy that we 
achieved. Parameters from which telescope delay 
variability could be determined were monitored 
during the observations (egs. ground weather con- 
ditions), and ionospheric data was available from 
GSFC (lata archives collected by Global Position- 
ing Satellites (GPS). 

The phase part of the complex visibility func- 
tion (called simply the phase) for each source a, 
can be written as 

<^f,„(t) = ^[B;,„.AK'^(t) + 

AB;,„(t) . K» + AC;,„(t) + ^Al^it)] + ^l^{t) (6) 

where 'Bi^m is the model separation between tele- 
scope I and m, AB(_m(t) is the separation er- 
ror; K" is the model position of the oth source. 



s www. sgl .crestech.ca/IVS-Analysis/ software-tools 
/calcsolve/dataiiles.htm 



A'K"'{t) is the position offset error for the ath 
source, ACi^m{t) is the residual clock delay be- 
tween telescopes, and AA'^^{t) is the residual tro- 
pospheric/ionospheric delay in the direction to the 
ath source. The structure of the source is given by 
the phase part of the complex spatial coherence 
function, ipfm{t). For this experiment, a has three 
values: a = (J0842-M835), a = 1 (J0839-M802), 
a = 2 (J0854-h2006). The occasional lobe am- 
biguities (arbitrary turns of phase) of the phase 
measurements were easily determined because of 
the accuracy of the delay model. 

3.2. The Source Structure 

The accuracy of the experiment depends in 
part on the stability of the three sources between 
September 4 to 12. In terms of variables in Eq. 
(6), the intrinsic change of position, AK."-{t), of 
each source, and the visibility structure phases 
variations, tj^f,-^(t), must be small. The struc- 
tures for each source on each day were determined 
from their respective complex visibility function 
using self-calibration techniques for radio imag- 
ing (Walker 1995). This imaging/deconvolution 
iteration process does not determine the accurate 
positions of each source, only its shape. The de- 
rived images of the three quasars during the ex- 
periment are shown in Fig. 1. All three sources 
show the typical structure associated with most 
luminous radio sources, a bright component (core) 
at one end of the structure containing an appre- 
ciable part of the emission, and more extended 
emission, often with a secondary bright compo- 
nent. There were no apparent changes in the in- 
tensity and shape of the sources, except for small 
changes in J0854-I-2006 at the 1% contour level. 
These properties are in agreement with previous 
observations of J0842 (sources names will be ab- 
breviated) and J0839 which show little long term 
structure changes (Fey & Chariot 2000; Beasley 
et al. 2001), whereas J0854 is a known variable 
source. 

The observation dates are symmetric with re- 
spect to September 8. Hence, even if there are 
long-term small structure changes over the exper- 
iment, the deviation of the average structure phase 
of each source, based on observations averaged on 
September 4,7,9 and 12, should be equal to first 
order to that on September 8 (this averaging is 
discussed below). Nevertheless, the change in the 



4 



measured phase associated with the structure dif- 
ference in J0854 between September 4 and 12 is 
less than 0.003 turn for any basehne, and this 
corresponds to an effective positional error of less 
than 3 /xarcsec, even without the above averaging 
which was incorporated in the reductions. An- 
other indication of very little change over the 7-day 
duration of the observations is that the separation 
of the radio core and the secondary component for 
all sources was stable to < 20 /iarcsec, an accuracy 
consistent with the sensitivity of the observations 
and the somewhat diffuse nature of the secondary 
peak. 

The large distance to the radio sources guar- 
antees that any proper motion will be much less 
than one microarcsec over the week period of the 
experiment. However, some compact components 
(< 50 /iarcsec) in radio sources vary slightly in 
intensity over hour time-scales because of galac- 
tic scattering at a distance of « 300 pc (Bignall 
et al. 2003). It is possible that such scattering 
also irregularly moves the apparent position of the 
component. However, we detected no short-term 
intensity variations in any source greater than 4% 
over periods of a few hours, and concluded that 
the compact components are not affected by in- 
terstellar scattering. However, we emphasize that 
such movement of the source position, if it existed, 
could be detected in this experiment as larger than 
expected residual phases with time. 

3.3. Removing the Temporal and Spatial 
PhEise Errors 

With the removal of the source structure terms 

(we used the average source structure over the five 
observing days), the phase in Eq. (6) becomes 
circular among the telescopes; that is <pi^rn = 
4>i.n + 4'n.rm apart from stochastic noise. There- 
fore, there are only 10 independent phases for an 
array of 11-telescopes. We can thus choose a ref- 
erence telescope, M, and write Eq. (6) for each 
source a and telescope Z, at time ta as 

^LMi*a)=fl,Mita) + 
gZ,M(ta,K'^)+AK'^(f,)-B;,M (7) 

where fi^uita) is a function of time, which is dom- 
inated by the clock drifts and other instrumental 
variations at each telescope, and is independent 
of the source position. The term g(,M(ia)K'*) de- 
scribes the temporal and spatial properties of the 



delay screen above each telescope near the posi- 
tion of the three sources. The last term in Eq. (7) 
cannot be written as part of the first two terms 
since it depends on an unknown, but small, offset 
for each source from the assumed value. We used 
Mauna Kea (MK), HI as the reference telescope 
for most of the observations since this telescope 
is at a high site with generally good atmospheric 
conditions. For the first three hours of each day 
when the source elevation at MK was < 20°, we 
used the Los Alamos, NM telescope as the refer- 
ence telescope. 

The properties of the several components of de- 
lay which contribute to the time variable phase 
screen over each telescope are well-known. The 
phase error associated with the uncertain location 
of a telescope AB; j^, which includes any devia- 
tion of the earth rotation and orientation and nu- 
tation terms from those used in correlator CALC 
Model, are smooth functions of time and direction 
in the sky. The component produced by the chang- 
ing refractivity of the troposphere and ionosphere 
toward the sources from relatively large clouds and 
diurnal changes persists for tens of minutes and 
are generally larger than 5° in angular size. On the 
other hand, small tropospheric clouds and iono- 
spheric events produce quasi-stochastic changes 
above each telescope. In their extreme behavior, 
these rapidly changing delays have time scales less 
than a few minutes and degrees, making astromet- 
ric VLBI observations useless. 

We can expand the first two terms in Eq. (7) 
into a Taylor series around (ta^^o) ^md (K"— K'^). 
Since J0842 and the two calibrators lie nearly lin- 
early in the sky, the angular part of the expan- 
sion is only needed in the direction of the source 
alignment in the sky*". The first order terms can 
be estimated by a sum of the observed calibrator 
phases, weighted by their separation from J0842 in 
the sky, and in time. Thus, the calibrated phase 
for J0842, $(to);,M at its observation time, to, for 
any day for telescope I becomes 

^{to)iM - 0°(io)/,M - {Q.m\to)i,M + Q.2Q4?{to)iM) (8) 
« B,,M-(AK°(io) - 0.80AKi(to) - 0.20AK2(fo)) 

where (/'"^(to) and (t)^{to) are the phases of sources 
1 or 2, each measured at times ti and t2, but 



Three calibrator sources are needed in the general case. 
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linearly-interpolated to time to- In other words, 
the algorithm in Eq. (8) defines how the measured 
phases for J0839 and J0854 are combined to re- 
move, to first order, the temporal and angular 
phase change in the sky at the position of J0842 at 
the time of its observation at to- The calibrated 
phase of J0842, thus, depends explicitly on the 
change of position of the three radio sources with 
time. The short time-scale and small angular-scale 
fluctuations arc not removed, but these arc nearly 
stochastic and average out over longer periods of 
time. They add to the uncertainty of the results, 
but do not contribute to any systematic biases. 
Second order terms, which are not removed from 
the interpolation used in Eq. (8), have been min- 
imized by using precise apriori parameters avail- 
able from the CALC package and available tropo- 
spheric and ionospheric data associated with each 
telescope from the GPS data base. 

Examples of the effectiveness of this calibration 
process are shown in Fig. 2 for the Owens Val- 
ley, CA to Mauna Kea, HI baseline on September 
9. The observed phases, (j>°'{t), (top plot) for the 
three sources follow each other over the day, but 
are separated at any time by an amount which is 
varying, caused by a phase wedge in the atmo- 
sphere. The displacement of phase between J0854 
and J0842 is larger and in the opposite sense than 
that from J0839 and J0842, as expected from lin- 
ear phase wedge. The calibrated phase $ of J0842, 
(bottom) is relatively smooth and near zero phase. 
The relatively small scatter is caused by fluctua- 
tions of phase shorter in temporal scale than about 
five minutes, or smaller in angular scale than five 
degrees. The signal to noise inherent in the obser- 
vations cause only a small part of the scatter. 

There are periods of time when the phase sta- 
bility shown by all three sources noticeably deteri- 
orates, and the calibrated phase of J0842 becomes 
large. These periods often occur when the source 
is at a low elevation, typically less than 20°, when 
large and non-linear tropospheric and ionospheric 
phase are more likely. Occasionally at high ele- 
vations during periods of rainy and windy condi- 
tions, VLBI observations are useless. These peri- 
ods of obvious poor phase stability were removed 
from further consideration, and consisted of 15% 
of the original data. On the average the phase sta- 
bility at sites with dry conditions was more stable 
than at those with humid conditions. 



The calibrated phases, $;.M(t), were then av- 
eraged for each day over a one hour period, with 
an estimated rms error determined from the scat- 
ter, which is about 0.02 turns per point. In Fig. 
3, we show these averaged calibrated phases for 
two telescopes. The calibrated phases repeat ex- 
tremely well from day to day to an rms scatter 
of about 0.02 turns. Even a slight negative off- 
set at GST» 15 hr repeats, and is produced by 
a '--^ 50 /xarcscc offset of the source position from 
that assumed in the model, as given in Eq. (8). 

3.4. The Deflection on September 8 

The difference between the calibrated phase on 
September 8, and the average of the phases from 
the other days, is shown in Fig. 4 for the same two 
telescopes as in Fig. 3. Assuming that the position 
of all three sources remain fixed with respect to 
the apriori model, which includes the radial and 
retarded gravitational bending predicted by GR, 
the phases in Fig. 4 should scatter around zero 
phase, corresponding to 5 = 0. The curve given 
by (5 = — 1 shows the expected calibrated phase on 
September 8 if the retarded deflection term were 
zero. Fig. 4 already demonstrates that (5 w 0, 
and that the experiment has clearly detected the 
retarded component of the light deflection. 

An analysis of the effect of the jovian magne- 
tosphere, has been described fully (Kopeikin & 
Fomalont 2002), and is summarized here. The 
Galileo spacecraft provided an estimate of the av- 
erage electron plasma density, A^'o, near the jovian 
surface. We assumed a spherically-shaped mag- 
netosphere with a plasma density which decreases 
as No{Rj/rf+^ (where Rj = 7.1 x lO"^ m is the 
mean radius of Jupiter and A is an exponent re- 
flecting different magnetospheric models), and in- 
tegrated the estimated plasma delay at 8.45 GHz 
along the propagation path of J0842 at closest ap- 
proach (3.7' or 13Rj). We obtain magnetospheric 
deflection estimates of 17.5, 1.0, 0.02 /iarcsec for 
A =0.0, 1.0, 2.0, respectively. Based on the anal- 
ogy with radial distribution of the solar corona, 
which suggests A « 0.33, we have determined the 
position shift associated with this jovian magne- 
tosphere model (it is in the opposite direction of 
the gravitational radial deflection of light), and we 
have shown its estimated phase contribution by 
the dashed line in Fig. 4. This estimate is prob- 
ably substantially larger than the actual value on 
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this day. 

There are two ways in which we can combine 
the data in order to determine the mean gravita- 
tional deflection of J0842. First, we averaged the 
data for each telescope, and determined 6 asso- 
ciated with each telescope, and these results are 
shown in Table 1, listed in decreasing order of ac- 
curacy. We find that the most accurate determina- 
tions are associated with telescopes in the south- 
west USA. about 5000 km from MK, because they 
are located at places with a relatively stable and 
dry atmosphere. Telescopes in more humid loca- 
tions give poorer results. The Effclsbcrg telescope, 
however, was critical in the determination of the 
accurate structure of the radio sources. All tele- 
scope determinations are consistent with 6 = 
when their estimated errors are considered. 

Another method of display is to produce an im- 
age from the calibrated phase data for all base- 
lines, which arc measured during any one-hour 
period. The location of the peak of the image 
is the position of J0842, and these positions are 
shown in Fig. 5. The east/west position of the 
source is more accurately determined than the 
north/south position since the array spans mostly 
in the east /west direction. The east /west posi- 
tions clearly show that the retarded component of 
the light deflection has been detected, and is con- 
sistent with GR. The red line shows the estimate 
of a reasonable jovian magnetospheric refraction, 
and is much smaller than the retarded deflection, 
particularly in the east/west direction. 

The weighted average value of the data for the 
observations over September 8, regardless of how 
the data are averaged (by telescope or by time), is 

6 = -0.02 ± 0.19 l-(7 error 

Hence, the measured retarded deflection associ- 
ated with the retarded position of Jupiter given by 
Eqs. (l)-(3), is in good agreement with the predic- 
tion of GR. Our interpretation of 6 sets a bound 
of Cg = (1.06 ± 0.21)c (see §4.b) 

4. Discussion 

4.1. The Measurement Results 

The measurement of the retarded deflection 
component in Eq. (1) is an accurate experimen- 
tal conflrmation of General Relativity in the solar 



system, using the motion of Jupiter. To deter- 
mine 6 with sufficient accuracy, improvements in 
the technique of phase-reference VLBI were devel- 
oped. 

The use of two calibrators, on either side of 
J0842, improved the positional accuracy of this 
phase- referencing VLBI experiment by a factor of 
more than three compared with previous obser- 
vations of this type. Because of the relative sim- 
plicity of the experimental concept as a classical 
deflection- type experiment, the uniform observa- 
tional procedure over the entire experiment, the 
data redundancy provided by the number of tele- 
scopes and the number of days, we are convinced 
that the result and error estimate are accurate and 
demonstrated from the Fig. 4 and Fig. 5 which dis- 
play a significant part of the data. The refraction 
from the jovian magnetosphere is at least as small 
as we have estimated (Kopeikin & Fomalont 2002). 
Any position jitter caused by interstellar scatter- 
ing is not significant since the scatter in Fig. 5 
around the expected retarded position is no larger 
than that expected from the errors derived from 
the scatter in the observations among all of the 
days. 

The two-calibrator technique, developed in this 
experiment to remove the delay variations of the 
tropospherically-induced phase wedge above the 
telescopes, can be used for determinations of the 
proper motion and parallax of radio sources. With 
the relative positional sensitivity of 10 /xarcsec, the 
trigonometric parallax of stars and pulsars as far 
as 10 kpc can be determined; spacecraft can be 
tracked with an accuracy of 30 m near Jupiter; 
and the proper motion of radio sources in M31 can 
be determined to 5 km/sec (1 /xarcsec/yr) over ten 
years. However, the significant structure changes 
in many calibrator sources over periods longer 
than a few months may dominate the error budget 
at angular levels of > 50 /^arcsec for the long-term 
astrometric projects. Expected improvements in 
the accuracy of the celestial frame over the next 
decade may permit the determination of calibra- 
tor core positions and changes needed to roach 
10 /xarcsec accuracy (Ma et al. 1998; Jacobs ct 
al. 2002) 

4.2. The Measurement Interpretation 

Our interpretation of this novel deflection ex- 
periment comes directly from the formulation of 
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Kopeikin (2001, 2003) where the retarded deflec- 
tion term (and the availabiUty of the September 8, 
2002 Jupiter experiment) was initially proposed. 
The interpretation is that the position of Jupiter 
at the VLBI delay given by Eq. (1) must be taken 
at the retarded time with respect to the time of ob- 
servation, as defined by the sohition of the gravity 
null cone Eqs. (2) and (3). This retardation is due 
to the finite speed of gravity and can be observed. 
The measiired retarded deflection, as expressed by 
the parameter S, is in agreement with GR to 20% 
accuracy, and that the formulation, summarized 
by Eqs. (l)-(3); is confirmed to the first order in, 
v/c. If (5 ^ 0, then the constant Cg in the wave 
operator of the Einstein equations is not equal to 
the speed of light numerically, and GR (the grav- 
ity null cone) would not be invariant with respect 
to the Lorentz transformation. 

Another interpretation has been proposed by 
Will (2003) who made use of a matching tech- 
nique that was first employed for solving the same 
problem by Klioner & Kopeikin (1992). This for- 
mulation assumes that the position of Jupiter in 
Eq. (1) must be taken at the time of the clos- 
est approach of light to Jupiter. Will's formula- 
tion is not Lorentz-invariant but approximates the 
Lorentz-invariant formulation of the time delay up 
to linearized terms of order v/c. 

Both formulations* give identical predictions 
for the deflection of light by a moving body to 
first order in (v/c). However, the physical mean- 
ing attached to the measurement is quite differ- 
ent. The Kopeikin (2001) interpretation is related 
to the Lorentz-invariant definition of the gravity 
null cone and the finite speed of propagation of 
gravity as predicted by the Lienard-Wiechert so- 
lution of Einstein's equations. Kopeikin (2003) 
also associates the v/c effects with the gravito- 
magnetic dragging of the light ray caused by the 
translational motion of Jupiter with respect to the 
barycenter of the solar system. The Will (2003) 
and Asada (2002) interpretation is that the v/c 
term is related to the aberration of light that 

'We suggest that the reader compare the treatment of the 
September 8, 2002 experiment by Kopeikin (2001, 2003), 
by Will (2003), by Asada (2002), for a fuller understand- 
ing of the different interpretations and mathematical tech- 
niques involved in derivation of the basic equations. Papers 
by Kopeikin & Schiifcr (1999); Klioner & Kopeikin (1992); 
Klioner (2003); Carlip (2000) may be valuable as well. 



diverges from our interpretation at terms of or- 
der v'^ /(?. Our measurement, unfortunately, is 
not sufficiently accurate to detect terms of order 
v"^ /c? which can distinguish between the existing 
interpretations experimentally. This could be a 
challenge for the next generation of space VLBI 
or other astrometric space missions like SIM or 
GAIA. Relativistic deflection terms of order v'^ /c^ 
for light grazing Jupiter's limb, may reach a mag- 
nitude of 1 /uarcsec and can be, in principle, ob- 
servable. 

What is clear is that if this experiment had 
measured 5 signiflcantly different from zero, then 
a problem would exist with the current formu- 
lation of General Relativity and, as such, any 
interpretation of the result would be ambiguous 
and depend on the theoretical generalization of 
the Einstein equations that was chosen. How- 
ever, this measurement of the deflection of light 
by Jupiter is in good agreement with General Rel- 
ativity and demonstrates experimentally that the 
Einstein Equations are Lorentz-invariant and the 
retarded position of Jupiter in Eq. (1) should be 
used in the determination of the deflection of light 
in accordance with Eqs. (2) and (3). 
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Table 1 
Telescope Solutions for S 



Telescope 


Baseline to MK 
(km) 


5 


Hancock, NH 


7500 


-0.02 ±0.25 


Los Alamos, NM 


4970 


-0.20 ±0.34 


Kitt Peak, AZ 


4490 


±0.52 ±0.38 


Owens Valley, CA 


4015 


-0.45 ±0.42 


Pic Town, NM 


4800 


+0.30 ±0.44 


Fort Davis, TX 


5130 


-0.09 ±0.52 


North Liberty, lA 


6160 


-0.37 ±0.53 


Brewster, WA 


4400 


+0.14 ±0.62 


Saint Croix, VI 


8610 


-0.82 ±0.85 


Effelsberg, Germany 


10300 


+1.94+1.60 


<S> 




-0.02 ±0.19 
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J0854+2006 

Fig. 1 - The Structure of the Three Quasars during the Experiment: The 

contour images for the three quasars on September 4, 8 and 12, are shown from left to 
right. North is up and cast is to the left. The first row shows J0842+1835 with a lowest 
contour of 3.0 mjy/beam, about 1% of the peak brightnesses of 306.4, 308.5, 308.4 
mjy/beam from left to right. The second row shows J0839+1802 with a lowest contour 
of 2.0 mJy/beam, about 2.8% of the peak brightnesses of 74.6, 75.5, 70.1 m.Jy/bcam. 
The third row shows .10854+2006 with a lowest contour of 7.0 m.ly, about 0.4% of 
the peak brightnesses of 1834, 1778, 1799 mJy/beam. The contour levels increase 
by a factor of two in brightness. The source J0854+2006 shows the emergence of 
a faint component to the east and a brightening of a weak component to the west. 
The resolution is 1.5 x 1.0 mas (full-width at half-maximum of an elliptical Gaussian 
beam), and the tick marks are separated by 1 mas (except 2 mas for the N/S axis of 
J0842+1835), and the lowest contour levels are about 5 times the rms noise. 
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Fig. 2 - Removing the Residual Tropospheric and Ionospheric Delay Vari- 
ations: The top plot shows the measured phase for the Owens Valley, CA, to the 
Mauna Kea, HI baseline on September 9 for each source. The red points are for 
J0842, the green points for J0839 and the purple points for J0854. The bottom plot 
shows the corrected phase for J0842 after using the linear interpolation of J0839 and 
J0854, defined by Eq. (7). The large phase gradient between UT 14h and 15h is pro- 
duced by the low elevation of the sources, and these data have been removed from 
further consideration. 




Fig. 3 - The Calibrated Phases for the HN-MK and the LA-MK baselines 
for all Days. The calibrated phases, averaged to one hour intervals on each day, for 
the Hancock, NH (top) and Los Alamos, NM (bottom) baselines to Mauna Kea, HI 
are connected by the line segments, with the day indicated at the right end of each 
line. The September 8 line is black. The error bars (not shown) are typically 0.02 
turns per point. The Hancock telescope was not in use on Sep 4 and 12 because of 
mechanical problems. 
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Fig. 4 - The Differential Calibrated Phase for September 8: The phase dif- 
ference between the calibrated phase on September 8, and the average of that on the 
other days are given for Hancock (top) and Los Alamos (bottom) telescopes to Mauna 
Kea. The errors are derived from the scatter among the off- Jupiter days and the esti- 
mated error per point on September 8. The blue curve for 5 = is that expected from 
GR; the pink curve for 5 = —1 is that expected for no retarded deflection component. 
The magnetosphere refraction estimate (see text) is shown by the dashed green line. 
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Fig. 5 - The Retarded Deflection of J0842+1835 on September 8, 2002: 

The plotted points and l-cr error estimates show the retarded position every hour for 
J0842+1835; that is, the measured deflection in the direction of Jupiter's motion. The 
radial deflection term of « 1200 /.tarcsec has been removed. The top plot shows the 
east /west position and the bottom plot shows the north/south position. The expected 
retarded deflection from GR {5 = 0) is shown by the thick purple line. The estimated 
jovian magnetosphere refraction is shown by the thin orange line. 



